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1.  INTRODUCTION 


This  interim  report  covers  the  period  June  15. 1990,  to  June  23,  1991 ,  and 
describes  the  scientific  progress  on  Contract  F1 9628-90-01 33  entitled  “Data 
Analysis  Suport  for  the  SPREE  Instruments  on  the  TSS-1  Shuttle  Flight."  The 
objective  of  this  contract  is  a  scientific  study  of  the  interactions  between  the  first 
flight  of  the  Tethered  Satellite  System  (TSS-1  )and  the  ionosphere  to  increase 
understanding  of  the  structure  and  current  characteristics  of  high-voltage 
sheaths  in  the  ionosphere. 

During  the  first  year  of  this  contract,  we  investigated  the  tether  current  and 
voltages,  the  shuttle  orbiter  ion  collection,  the  ion  current  detected  by  the 
SPREE  instrument,  the  effect  of  the  subsateliite  thrusters  on  the  neutral 
densities,  and  provided  support  for  the  TSS-1  program. 

The  technical  staff  who  contributed  to  the  research  described  in  this  report  are  I. 
Katz,  J.  R.  Lilley,  Jr,  T.  T.  Luu,  and  V.  A.  Davis. 
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2.  CALCULATION  OF  TSS-1  USING  EPSAT  (SPECIAL  VERSION) 


2.1  INTRODUCTION 

Calculations  of  the  tether  current  and  voltages  over  the  mission  were  performed 
using  a  specially  modified  version  of  EPSAT  (Environment-Power  System 
Analysis  Tool).  The  model  consisted  of  an  arbiter  with  the  subsatellite  extended 
20  km  out  of  the  bay.  Collection  of  electrons  was  the  product  of  the  1  -sided 
thermal  flux  and  2nr^  times  the  minimum  of  the  space-charge-limited  or 
magnetic-limited  collection  radii.  The  ion  collection  is  just  the  space-charge- 
limited  ion  collection  to  a  sphere  with  effective  area  that  of  the  engine  bells. 

The  first  calculations  examined  the  effect  of  subsatellite  paint  resistance  on  the 
tether  currents.  It  is  uncertain  whether  the  subsatellite  paint  will  be  8K  Q,  the 
highest  impedance  in  the  tether  circuit,  or  whether  it  will  have  negligible  bulk 
resistance  at  the  operating  potentials.  With  the  core  gun  on,  using  a  perveance 
of  3.8  X  10'®,  tether  currents  were  calculated  for  the  mission  duration  of  different 
sphere  and  tether  resistances.  The  results  were  surprisingly  insensitive  to  the 
paint  resistance.  When  the  total  resistance,  tether  plus  subsatellite  paint,  was 
10K  fi,  the  peak  current  was  about  60  percent  of  that  for  an  assumed  2K  Q 
resistance.  Peak  subsatellite  potentials  differed  little,  because  they  occur  at  low 
plasma  density  when  low  ionospheric  currents  cause  the  resistive  voltage  drop 
to  be  negligible. 

The  other  sequence  of  EPSAT  studies  examined  orbiter  potentials  when  the 
tether  is  connected  to  the  orbiter  chassis  through  a  large  resistance.  The 
resistances  studied  were  between  25K  Q  and  350K  Q.  The  calculations  showed 
that  the  orbiter  would  charge  only  at  low  plasma  density,  again  because  the 
resistive  voltage  drop  is  small  due  to  low  tether  currents.  However,  orbiter 
potentials  as  large  as  -2200  V  resulted  when  the  total  tether,  subsatellite,  and 
shunt  resistance  was  35K  Q. 
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Figure  2.1  The  EPSAT  model  of  the  shuttle  orbiter. 
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Figure  2.2  Side  view  of  the  EPSAT  Orbiter  model. 


The  model  was  oriented  with  the  subsatellite  away  from  the  earth. 
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Max  Tother  Currant  vs.  Total  Reelstance 


Figure  2.5 


Results  for  a  low  current  time  on  orbit. 
No  gun,  35K  Q  impedance. 
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3.8E-6  Perveance,  10K  Q  Impedance. 
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Plasma  conditions  at  one  time  during  the  mission. 
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Figure  2.6 
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Potentials  With  35K  Ohm  Shunt 


Figure  2.7 

Orbiter  Potential  va.  Shunt  Resistance 


3.  CALCULATIONS  OF  ION  CURRENT  COLLECTED 
BY  SHUTTLE  ORBITER 


Using  NASCAP/LEO,  the  collected  ion  (oxygen)  currents  were  calculated  as  a 
function  of  voltage,  plasma  density,  and  the  orientation  of  the  shuttle  orbiter  with 
respect  to  the  ram.  The  ion  currents  were  found  by  varying  the  conducting 
surface  potentials  (the  engine  nozzles),  computing  the  ion  wake  formed  behind 
the  moving  shuttle,  calculating  the  resulting  spatial  plasma  potentials,  and  then 
pushing  ions  from  an  equipotential  plasma  sheath. 

The  shuttle  model  was  made  using  PATRAN.  The  model  coordinate  system  is  such 
that  the  tail  rises  in  the  +x  direction,  the  nose  is  in  the  +z  direction,  and  starboard 
side  of  the  shuttle  is  the  +y  direction.  For  the  purposes  of  this  calculation,  only  the 
engines  nozzles  were  considered  conducting  (aluminum  in  the  figure).  The 
calculations  used  a  nested,  cubical  grid  with  a  mesh  size  of  1 .645  meters.  The 
plasma  temperature  for  all  of  the  calculations  was  0.1  eV.  The  surface  potentials  of 
the  insulating  surfaces  were  defined  to  be  -0.5  volts  to  simplify  the  calculations. 

The  varied  parameters  were  the  potential  of  the  conducting  surfaces  (-10,  -100, 
-300,  and  -1000  volts),  the  piasma  density  (10**11  and  10**12  ions/meter**3), 
and  the  shuttle  orientation  with  respect  to  its  velocity  7500  meters/sec  (cargo 
bay  45  degrees  into  the  ram  and  into  the  wake).  The  corresponding  shuttle 
velocity  vectors  in  model  coordinates  are  (5303,  0,  -5303)  meters/sec  for  the 
bay  facing  in  the  ram  direction  and  (-5303,  0,  5303)  meters/sec  for  the  bay 
facing  in  the  wake  direction. 

The  results  of  the  calculations  follow: 


Bay  into  the  ram  (shuttle  velocity  was  (5303,  0,  -5303)  meters/sec) 
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Bay  into  the  wake  (shuttle  velocity  was  (-5303,  0,  5303)  meters/sec) 
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4.  CALCULATIONS  OF  ION  CURRENT  COLLECTED 
BY  THE  SPREE  INSTRUMENT 


Using  NASCAP/LEO,  ion  (oxygen)  currents  were  calculated  as  a  function  of 
voltage  and  the  orientation  of  the  SPREE  instrument  with  respect  to  the  flowing 
plasma.  The  currents  calculated  were  the  total  current  to  each  of  the  detector 
entrance  screens,  the  current  that  entered  the  detector  within  10  degrees  of  the 
surface  normal,  and  the  weighted  (by  particle  current)  average  angle  of 
particles  striking  the  detector  entrance. 

The  ion  currents  were  found  by  varying  the  surface  potentials  of  the  detector 
cases,  computing  the  ion  wake  formed  by  the  moving  instrument,  calculating  the 
resulting  spatial  plasma  potentials,  and  then  pushing  ions  from  an  equipotential 
plasma  sheath.  The  particle  velocity  and  current  were  used,  saved,  and 
postprocessed  to  compute  the  calculated  currents. 

The  SPREE  model  was  made  using  PATRAN.  The  model  coordinate  system  is 
such  that  the  detectors  rise  in  the  +y  direction  with  the  base  plate  in  the  x-z 
plane.  The  detector  with  screen  2  is  higher  in  z  than  the  detector  with  screen  1 . 
The  potentials  were  applied  to  the  case  of  the  detectors  and  their  screens.  The 
calculations  used  a  nested,  cubical  grid  with  a  mesh  size  of  0.04  meters.  The 
plasma  temperature  for  all  of  the  calculations  was  0.1  eV.  The  surface  potentials 
of  all  of  the  surfaces  other  than  the  case  of  the  detector  were  defined  to  be  -0.5 
volts  to  simplify  the  calculations. 

The  varied  parameters  were  the  potential  of  the  conducting  surfaces  (-100, 

-300,  and  -1000  volts)  and  the  SPREE  orientation  with  respect  to  the  flowing 
plasma  (or  orbiter  velocity  of  7500  meters/sec).  Two  orientations  were  studied. 
In  case  1 ,  the  flow  direction  corresponds  to  the  orbiter  nose  up  45  degrees  and 
the  cargo  bay  into  the  ram.  Case  2  arises  when  the  orbiter  is  oriented  with  the 
wings  parallel  to  the  flow  (shuttle  moving  to  port  or  starboard).  The 
corresponding  SPREE  velocity  vectors  in  model  coordinates  are  (0,  5303, 

-5303)  meters/sec  for  case  1  and  (7500,  0,  0)  meters/sec  for  case  2. 


The  results  of  the  calculations  follow: 


Case  1  Nose  up,  cargo  bay  into  the  ram.  (SPREE  velocity  was  (0,  5303,  -5303) 
meters/sec.)  In  this  orientation,  the  low-Z  detector  is  the  aft  detector. 


AFT  DETECTOR 


Voltage  (volts) 

Total  current  (nA) 

-100 

4.88 

1.15 

•1000 

1.61 

Current  within 

1 0  degrees  (^lA) 

.0303 

.055 

0. 

Weighted  average 
angle  detected  (degrees) 

43.8 

64.9 

59.4 

FORWARD 

DETECTOR 

Voltage  (volts) 

Total  current  (jiA) 

-JLflil 

2.43 

ilAi} 

.694 

-1000 

2.71 

Current  within 

1 0  degrees  (m-A) 

.0962 

0. 

.073 

Weighted  average 
angle  detected  (degrees) 

43.4 

49.3 

41.7 

Case  2  Orbiter  moving  sideways.  (SPREE  velocity  was  (7500,  0,  0) 

meters/sec.)  In  this  orientation,  the  low-Z  detector  is  looking  towards 
its  wake  and  the  high-Z  detector  is  looking  upstream. 


LOW-Z  DETECTOR  (DETECTOR  1),  WAKE  VIEW 


Voltage  (volts) 

-100 

Total  current  (nA) 

5.14 

4.93 

3.84 

Current  within 

10  degrees  (^lA) 

.102 

0. 

0. 

Weighted  average 

angle  detected  (degrees) 

47.5 

67.0 

56.7 

HIGH-2  DETECTOR 

(DETECTOR  2), 

RAM  VIEW 

Voltage  (volts) 

lUifl 

-1000 

Total  current  (nA) 

2.43 

5.18 

4.03 

Current  within 

10  degrees  ()iA) 

.098 

0. 

0. 

Weighted  average 

angle  detected  (degrees) 

43.1 

72.6 

59.4 
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5.  OBSTRUCTION  OF  CURRENT  OF  SPREE  INSTRUMENTS 
DUE  TO  RMS  AND  CAMERA 


5.1  OBJECTIVES 

This  chapter  describes  the  investigation  of  the  degree  of  obstruction  to  the 
SPREE  instrument  due  to  the  RMS  (Remote  Manipulation  System)  and  the 
camera  mounted  on  top  of  the  RMS. 


5.2  PROCEDURE 

The  calculations  were  done  using  NASCAP/LEO,  which  has  been  modified  to 
treat  wake  and  ram  particles  more  efficiently.  The  currents  calculated  were 

1 .  Total  current  to  the  SPREE  instruments  (the  two  boxes); 

2.  Total  current  to  each  of  the  detector  entrance  screens; 

3.  Current  that  entered  the  detector  within  1 0°  of  the  surface  normal;  and 

4.  The  weighted  (by  particle  current)  average  angle  of  particles  striking  the 
detector  entrance. 

5.3  RUN  CASES 

Ion  currents  to  the  detectors  were  calculated  for  three  cases: 

1.  SPREE  only  (repeat  previous  calculation); 

2.  SPREE  and  RMS  and  nonconducting  camera; 

3.  SPREE  and  RMS  and  conducting  camera. 
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5.4  PARAMETERS 


Parameters  common  to  all  three  cases  are  as  follows: 


Plasma  temperature 

0.1 

eV 

Ion  density 

10” 

m-3 

Orbiter  velocity 

7500 

m/s 

Inner  grid  mesh  size 

0.06 

m 

In  all  cases,  a  potential  of  -300  V  is  applied  to  the  surfaces  of  the  detector  cases. 
In  case  3,  the  same  potential  is  applied  to  the  camera  surfaces  as  well. 


5.5  ORIENTATIONS 


In  all  cases,  the  low-Z  detector  is  looking  toward  its  wake  and  the  high-Z 
detector  is  looking  upstream.  Different  orientations  used  in  the  calculations  are 
shown  in  table  5.1. 

Table  5.1  Detector  orientations  used  in  calculations. 


Case 

VelocitY(m/s) 

1 

0 

-10 

7386.0.-1302 

2 

10 

-10 

7274.1302.-1283 

3 

20 

-10 

6941.2565.-1224 

4 

0 

0 

7500.0.0 

5 

10 

0 

7386.1302.0 

6 

20 

0 

7048.2565.0 

7 

0 

10 

7386,0.1302 

8 

10 

10 

7274,1302,1283 

9 

20 

1.0 

6941,2565,1224 

The  symbols  0  and  (Jjare  defined  in  figure  5.1. 
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Figure  5.1  Definition  of  6  and  <|>. 


87 


88 


igure  5.2 


5.6  RESULTS 


All  currents  are  in  amperes.  The  low-Z  detector  is  facing  the  wake,  and  the 
high-Z  detector  is  facing  the  ram. 


Table  5.2  Total  currents  collected  by  both  detector  boxes. 


Orientation 

e 

0 

Case  1 

Case  2 

Case  3 

1 

0 

-10 

2.27e-4 

l.75e-4 

l.76e-4 

2 

10 

-10 

2.14e-4 

l.6le-4 

1.63e-4 

3 

20 

-10 

2.03e-4 

l.7le-4 

1.30e-4 

4 

0 

IQIII 

2.25e-4 

1.85e-4 

1.86e-4 

5 

10 

0 

2.03e-4 

1.52e-4 

1.51e-4 

6 

20 

0 

2.02e-4 

1.72e-4 

l.31e-4 

7 

0 

10 

2.25e-4 

1.83e-4 

1.83e-4 

8 

10 

10 

2.13e-4 

1.61e-4 

1.58e-4 

9 

■B 

10 

2.12e-4 

i.7le-4 

1.38e-4 

Table  5.3  Total  currents  collected  at  the  detector  screens. 


Low-Z 

Hiqh-Z 

Orientation 

6 

n 

Case  1 

Case  2 

CaSG  3 

Case  1 

Case  2 

Case  3 

1 

0 

-10 

6.82e-6 

7.31e-6 

9.52e-6 

8.15e-6 

7.72e-6 

6.86e-6 

2 

10 

-10 

6.43e-6 

6.39e-6 

6.79e-6 

5.84e-6 

4.77e-6 

6.55e-6 

3 

20 

-10 

5.69e-6 

6.63e-6 

4.76e-6 

5.76e-8 

3.90e-6 

5.60e-6 

4 

n 

0 

7.17e-6 

6.52e-6 

8.41e-6 

6.06e-6 

8.20e-6 

7.59e-6 

5 

10 

n 

5.41e-6 

5.79e-6 

6.17e-6 

6.05e-6 

6 

20 

■i 

5.16e-6 

6.28e-6 

4.17e-6 

5.25e-6 

5.45e-6 

6.53e-6 

7 

wm 

10 

6.06e-6 

6.13e-6 

6.67e-6 

6.32e-6 

6.69e-6 

6.29e-6 

8 

10 

10 

5.01e-6 

2.69e-6 

5.68e-6 

5.94e-6 

4.95e-6 

9 

20 

10 

4.52e-6 

5.l1e-6 

4.26e-6 

6.09e-6 

6.39e-6 

3.49e-6 

89 


Table  5.4  Total  currents  entering  the  detectors  within  10^ 
of  the  surface  normal. 


Low-Z 

Hiah-Z 

Orientation 

El 

n 

Case  1 

Case  2 

Case  3 

Case  1 

Case  2 

Case  3 

1 

0 

-10 

1.45e-7 

1.42e-8 

1.64e-7 

5  27e-7 

0. 

2 

10 

-10 

1.37e-7 

1.08e-8 

3.97e-8 

5  47e-7 

7.06e-7 

0. 

3 

-10 

1.34e-8 

4.04e-9 

6.40e-9 

6.91e-7 

7.87e-7 

1.73e-7 

4 

0 

3.94e-8 

2.05e-8 

4.53e-8 

2.61e-7 

9.23e-7 

5 

10 

0 

1.77e-8 

3.05e-10 

4.92e-8 

1  77e-8 

1.63e-7 

8.26e-7 

6 

20 

0 

4.16e-9 

1.09e-7 

3.02e-7 

7 

0 

m 

4.15e-8 

3.37e-10 

3.65e-9 

0. 

3.73e-7 

8 

10 

■1 

2.13e-20 

3.87e-20 

7.53e-8 

0. 

9.67e-8 

9 

20 

m 

1.37e-9 

1.34e-10 

2.47e-8 

5.65e-8 

0. 

6.89e-10 

Table  5.5  Weighted  average  angle  of  the  particles  striking  the  detectors. 


Low-Z 


Case  2 


50.8 


50,0 


48.0 


50.3 


48.1 


48.1 


Cases 


38.9 


41.4 


45.0 


42.3 


42.3 


Hiqh-Z 

Case  1 

Case  2 

Cases 

34.6 

36.9 

44.5 

27.8 

31.0 

23.8 

23.2 

38.0 

29.1 

33.4 

50.7 

26.5 

35.2 

34.0 

25.9 

31.9 

35.2 

28.0 

37.7 

51.9 

25,2 

37.1 

48.3 

23.8 

52.4 

8 


10 

20 


47.0 

54.9 

48.3 


47.4 


40.1 

40,7 


SPREE.  VeU(I.O.O).  rho=1.e11 


Figure  5,3  Sheath  for  case  1 — SPREE  only. 
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Figure  5.4  ParJcle  trajectories  for  case  1 — SPREE  only — projected  on 
the  (a)  X-Y  plane,  (b)  Y-Z  plane,  and  (c)  X-Z  plane. 
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Figure  5. 5  Sheath  for  case  2 — SPREE  and  RMS  and  nonconducting  camera. 


1 7.00< _  SPREE,  high  Z  detector,  Vel. =(1,0,0),  rho=1e11 

16.00^  Tu  II///  /X 


Figure  5.6  Particle  trajectories  for  case  2— SPREE  and  RMS  and  nonconducting 
camera  projected  on  the  (a)  X-Y  plane,  (b)  Y-Z  plane,  and  (c)  X-Z  plane. 
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Figure  5.7  Sheath  for  case  3 — SPREE  and  RMS  and  conducting  camera,  (a)  Cutplane 
at  X  =  13,  (b)  cutplane  at  Y  =  12,  and  (c)  cutplane  at  Z  =  10.5. 
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6.  CALCULATIONS  OF  ION  CURRENT  COLLECTED 
BY  SHUTTLE  ORBITER 


This  attachment  describes  the  calculations  of  the  collected  currents  as  a 
function  of  the  orientation  of  the  orbiter  with  respect  to  the  ram,  the  potential, 
and  the  plasma  density. 


6.1  ION  CURRENT  COLLECTED  AS  A  FUNCTION  OF  ORIENTATION 

The  calculations  were  done  using  NASCAP/LEO,  which  has  been  modified  to 
treat  wake  and  ram  particles  more  efficiently. 

The  following  procedure  was  used  in  this  set  of  calculations; 

1 .  Select  an  orientation  of  the  orbiter; 

2.  Compute  the  ion  wake  formed  behind  the  moving  shuttle: 

3.  Calculate  the  resulting  spatial  plasma  potentials; 

4.  Push  ions  from  an  equipotential  plasma  sheath  to  compute  the  collected 
current. 

This  set  of  calculations  was  done  with  the  following  parameters: 


Plasma  temperature 

0.1 

eV 

Ion  density 

10^2 

m-3 

Orbiter  velocity 

500 

m/s 

Conductor  potential 

-300 

V 
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Figure  6.1  Shuttle  orientations  used  in  these  calculations. 


As  shown  in  figure  6.1 ,  the  nose  of  the  orbiter  was  tilted  45°  up,  and  the  arbiter 
was  rotated  around  the  vertical  axis.  The  arbiter's  belly  was  facing  the  ram  at  0° 
rotation  and  the  wake  at  1 80°.  The  rotation  step  was  1 5°.  Shuttle  velocity 
vectors  (traveling  at  7500  m/s)  are  given  in  table  6.1 : 


Table  6.1  Velocity  vectors  of  the  shuttle  in  the  different  orientations. 


Case 

Anqlen 

Velocitv(m/s) 

1 

0 

-5303.  0.  5303. 

2 

15 

-5123  1941.  5123. 

3 

30 

-4593.  3750.  4593 

4 

45 

-3750.  5303.  3750. 

5 

60 

-2652.  6495.  2652. 

6 

75 

-1373.  7244.  1373. 

7 

90 

0.  7500.  0. 

hbb 

1373.  7244.  -1373. 

HIS 

2652.  6495.  -2652. 

10 

135 

3750  5303.  -3750. 

1 1 

150 

4593.  3750.  -4593. 

12 

165 

5123  1941.  -5123. 

13 

180 

5303.  0.  -5303. 
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Figure  6.2  shows  the  model  of  the  shuttle  used  for  these  calculations. 


(a) 

Figure  6.2  Patran  model  of  the  shuttle  from  three  directions. 
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6.2  RESULTS 


At  0°  angle,  the  orbiter  is  traveling  with  its  belly  facing  the  ram  directly  and  the 
engine  nozzles  were  in  the  wake.  Thus,  little  current  was  collected. 

As  the  orbiter  rotates,  more  of  the  engine  nozzles  were  exposed  to  the  ram  side 
of  the  orbiter  and  more  current  was  collected. 

At  90°,  the  orbiter  is  moving  sideways  to  the  ram,  and  most  of  the  sides  of  the 
nozzles  were  exposed  to  the  ram.  The  currents  to  the  nozzles'  sides  peaked  at 
this  point  and  remained  steady. 

At  180°,  the  nozzles  faced  the  ram  direction,  and  the  currents  to  the  faces  of  the 
engine  nozzles  peaked. 

Table  6.2  and  figure  6.3  show  ion  current  at  different  rotation  angles.  All 
currents  are  in  amperes. 


Table  6.2  Collected  current  at  different  orientations. 


Case 

Angle 

Exterior 

Nozzle 

Face 

Bay 

Inside 

Doors 

Nozzle 

Doors 

Total 

1 

1.1e-6 

4.2e-3 

O.Oe+0 

O.Oe+0 

5.5e-3 

9.675e-3 

2 

15 

1.1e-5 

4.3e-3 

O.Oe+0 

O.Oe+0 

4.7e-3 

9.086e-3 

3 

30 

2.0e-5 

4.3e-3 

O.Oe+0 

O.Oe+0 

5.7e-3 

1.00le-2 

4 

45 

1 .8e-5 

4.4e-3 

O.Oe+0 

O.Oe+0 

9.2e-3 

1 .360e-2 

5 

60 

1.1e-5 

6.2e-3 

O.Oe+0 

O.Oe+0 

l.7e-2 

2.347e-2 

6 

75 

5.4e-5 

8.0e-3 

O.Oe+0 

O.Oe+0 

2.6e-2 

3.418e-2 

7 

90 

1 .4e-4 

l.0e-2 

O.Oe+0 

O.Oe-t-O 

3.0e-2 

4.082e-2 

8 

105 

3.4e-4 

1.4e-2 

O.Oe+0 

O.Oe+0 

3.2e-2 

4.655e-2 

9 

120 

6.5e-4 

1 .9e-2 

O.Oe+0 

O.Oe+0 

3.3e-2 

5.247e-2 

10 

135 

5.9e-4 

2.3e-2 

O.Oe+0 

O.Oe+0 

3.3e-2 

5.71 1e-2 

11 

150 

6.3e-4 

2.7e-2 

O.Oe+0 

O.Oe+0 

3.3e-2 

6.071  e-2 

12 

165 

5.1e-4 

3.0e-2 

O.Oe+0 

O.Oe+0 

3.2e-2 

6.313e-2 

13 

180 

2.7e-4 

3.le-2 

O.Oe+0 

O.Oe+0 

3.3e-2 

6.406e-2 
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Figure  6.3  Ion  currents  vs.  rotation  angles. 

6.3  ION  CURRENT  COLLECTED  AS  A  FUNCTION 
OF  POTENTIAL  AND  PLASMA  DENSITY 

The  following  procedure  was  used  in  this  set  of  calculations: 

1 .  Compute  the  ion  wake  formed  behind  the  moving  shuttle. 

2.  Calculate  the  resulting  spatial  plasma  potentials. 

3.  Push  ions  from  an  equipotential  plasma  sheath  to  compute  the  collected 
current. 

The  varied  parameters  are  conductor  potentials,  orientation,  and  plasma 
density. 

1 .  Conducting  potentials:  -1 0,  -1 00,  -300,  -1 000  V 

2.  Velocity  vectors:  (5303,0,-5303)  for  bay-to-ram  direction  and 
(-5303,0,5303)  for  belly-to-ram  direction. 
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Tables  6.3  through  6.6  give  the  results  of  the  calculations.  All  currents  arc  ,r. 
amperes. 


Table  6.3  Ion  current  collected  with  bay-to-ram  and  10’ Vm^  plasma  density. 


Pot. 

TefI 

Silv 

Kapt 

Si02  i 

Alum  ^  Total 

10 

1.2e-04 

2.1e-03 

O.Oe+00 

O.Oe+00  i 

1.5e-03 

3.721e-03 

-100 

7.9e-05 

3.3e-03 

O.Oe+00 

4.1e-03 

7.414e-03 

-300 

6.5e-05 

■HRiltCTi 

O.Oe+OO 

O.Oe+OO 

6.1e-03 

1.011e-02 

-1000 

1.0e-04 

4.7e-03 

O.Oe+00 

O.Oe+00 

9.9e-03 

1.466e-02 

Table  6.4  Ion  current  collected  with  bay-to-ram  and  10’ ^/m^  plasma  density. 


Pot.  : 

TefI 

Silv 

Kapt 

Si02 

Alum 

Total 

1 

o 

1.2e-03 

2.1e-02 

O.Oe+OO  * 

O.Oe+OO 

1.5e-02 

3.703e-02 

-100 

6.9e-04 

2.5e-02 

O.Oe+OO 

2.1e-02 

4.618e-02 

-300 

3.1e-02 

O.Oe+OO  ! 

O.Oe+OO 

3.3e-02 

-1000 

4.8e-04 

3.8e-02 

O.Oe+OO  i 

O.Oe+OO 

5.1e-02 

8.916e-02 

Table  6.5  Ion  current  collected  with  beliy-to-ram  and  10”/m3  plasma  density. 


Pot. 

TefI 

Silv 

Kapt 

Si02 

Alum 

Total 

-10 

1.5e-06 

5.0e-07 

2.2e-07 

2.5e-05 

2.877e-05 

-100 

3.6e-07 

0  Oe+00 

O.Oe+OO 

9.3e-04 

1.794e-03 

-300 

6.4e-05 

2.ie-03 

O.Oe+OO 

O.Oe+OO 

3.992e-03 

-1000 

4.3e-05 

le-03 

O.Oe+OO 

O.Oe+OO 

4.0e-03 

9.160e-03 

Table  6.6  Ion  current  collected  with  belly-to-ram  and  10’ ^/m^  plasma  density. 


Pot. 

TefI 

Silv 

Kapt 

Si02 

Alum 

Total 

-10 

2.7e-05 

1.8e-05 

1.0e-07 

6.0e-06 

2  5e-04 

2.992e-04 

-100 

i.8e-09 

3.3e-08 

7.6e-04 

9.6l7e-04 

-300 

1. IE-06 

4.2E-03 

O.OE+OO 

5.5E-03 

9.675E-03 

-1000 

3.5E-04 

1  5E-02 

O.OE+OO 

3.2E-08 

1.6E-02 

3.147E-02 

1  1  0 


7.  NEUTRAL  DENSITIES  ABOUT  TSS-1  SUBSATELLITE 


We  used  EPSAT  to  examine  the  plasma  environment  and  sheath  structure  of  the 
subsatellite  when  the  subsatellite  is  at  the  highest  potential  it  reaches  with  the 
electron  gun  on.  The  perveance  of  the  gun  is  taken  to  be  3.8  x  10“®,  and  the  net 
resistance  of  the  subsatellite  paint,  tether  wire,  and  shunt  is  tal  en  to  be  10  MQ. 

Figure  7.1  shows  the  subsatellite  potential  as  a  function  of  time.  The  peak 
potential  is  3610  V  at  1 125  s  into  the  mission.  The  plasma  environment  at  this 
time  is  shown  in  figure  7.2,  and  the  potentials  are  shown  in  figure  7.3.  Table  7.1 
shows  the  variation  in  potential  with  time  near  the  peak  potential  time.  Table  7.2 
shows  the  Parker-Murphy  (magnetic-limited)  and  space-charge-limited  radii  and 
collected  currents  at  the  peak  potential. 

Subsateilite  Potentlci  (Gur  On) 


time  into  mission  (sec) 


Figure  7.1  The  subsatellite  potential  as  a  function  of  time.  The  perveance  of 
the  gun  is  3.8  x  10“®,  and  the  net  resistance  of  the  subsatellite 
paint,  tether  wire,  and  shunt  is  10  Mt2. 
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Figure  7.2  EPSAT  screen  that  shows  the  plasma  environment  at  the  high 
potential  point  in  the  orbit. 
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Figure  7.3  EPSAT  floating  potential  screen.  This  screen  summarizes  the 
environment,  the  potentials  on  each  object,  and  the  collected 
current  for  each  object. 
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Table  7.1  Variation  of  floating  potential  with  time. 


Time  (s) 

Subsatellite  Potential  (V) 

1100 

3607 

1123 

3610 

1200 

3549 

1300 

3271 

1400 

2779 

Table  7.2  Current  collected  by  the  subsatellite  at  3609  V 
in  the  above  environment. 

Collection  Radius  (m)  Current  Collected  (A) 

Parker-Murphy 

(Magnetic-Limited)  2.88  -0.0305 

Space-Charge-Limited  6.13  -0.276 


We  also  used  EPSAT  to  examine  the  neutral  densities  and  the  likelihood  of  bulk 
ionization  under  conditions  of  high  outgassing,  the  2  N  in-line  thruster,  and  the 
2.5  kg/hr  leak. 

In  order  to  determine  if  enough  neutral  atoms  will  be  outgassed  from  the  subsatellite 
to  present  a  risk  of  bulk  ionization,  we  examined  outgassing  densities  for  the  case 
where  all  the  surfaces  emit  0.5  W  m“2.  This  is  the  outgassing  rate  of  G-10  after 
exposure  to  air  and  provides  an  upper  bound.  Figures  7.4(a)  and  (b)  show  the 
neutral  densities  about  the  subsatellite  due  to  outgassing  at  this  high  rate. 

Figure  7.5  is  the  neutral  column  density  EPSAT  screen.  The  screen  shows  that  the 
column  density  of  outgassed  neutrals  from  a  point  just  outside  the  subsatellite  (a 
1.02  m  cube)  to  a  point  2.88  m  away  (the  Parker-Murphy  radius)  is  1.6  x  1017.  vvith 
a  cross  section  of  3  x  1 0~20^  the  ionization  fraction  (number  of  ionizations  as  a 
fraction  of  the  number  needed  to  cause  ionization)  is  0.0048.  EPSAT  (ionization 
screen)  estimates  the  ionization  fraction  to  be  0.15  when  the  space-charge  sheath 
size  is  used  as  the  average  path  length.  Since  breakdown  is  expected  when  the 
ionization  fraction  approaches  1  and  the  outgassing  rate  used  is  orders  of 
magnitude  higher  than  expected,  no  bulk  ionization  is  expected  from  outgassing 
alone. 
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Figure  7.4  Neutral  densities  about  the  subsatellite  from  outgassing  at 

0.5  W  m-2,  which  is  a  maximal  rate,  (a)  is  a  3-dimensional  plot, 
and  (b)  is  a  contour  plot.  The  units  of  position  are  meters  and  of 
density  are  molecules  per  cubic  meter. 
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Figure  7.5  EPSAT  neutral  column  density  screen  for  outgassing  only. 

We  used  the  EPSAT  nozzle  model  to  model  one  nozzle  of  the  2  N  in-line  thmster  on 
the  subsatellite.  Figure  7.6  shows  the  parameters  used  in  the  model.  Figures  7.7 
show  the  gas  density  about  the  subsatellite  with  operation  of  one  nozzle  of  the  2  N 
in-line  thruster  only.  We  computed  column  densities  for  the  thruster  gas.  The 
computations  are  shown  in  figure  7.8.  The  column  density  for  a  2.88  m  column  offset 
from  the  axis  of  the  nozzle  by  0.51  m  is  5  x  10^0.  The  column  density  for  a  2.88  m 
column  normal  to  the  axis  of  the  nozzle  starting  0.51  m  from  the  nozzle  is  6  x  lO^^. 
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Figure  7.6  Nozzle  model  screen  showing  parameters  used  to  model  one 
nozzle  of  the  2  N  in-line  thruster. 
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7.7  Gas  densities  due  to  one  nozzle  of  the  2  N  in-line  thruster,  (a)  is  a  contour  plot,  (b)  and  (c)  show 
the  gas  density  along  the  thruster  axis  and  normal  to  the  nozzle  axis.  All  distances  are  in  meters 
and  all  gas  densities  are  in  molecules  per  cubic  meter. 
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Figure  7.8  Column  density  screens  for  the  computation  of  the  column  density 
(a)  along  the  nozzle  axis  and  (b)  normal  to  the  nozzle  axis. 

Nozzle  parameters  for  the  other  thruster  nozzles  were  also  determined.  EPSAT 
screens  describing  each  of  the  nozzles  are  shown  in  figure  7.9  through  7.12. 
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Figure  7  9  EPSAT  screen  showing  nozzle  parameters  for  one  nozzle  of  the  1 
N  in-line  thruster. 
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Figure  7.10  EPSAT  screen  showing  nozzle  parameters  for  one  nozzle  of  the 
in-plane  thruster. 
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Figure  7.11  EPSAT  screen  showing  nozzle  parameters  for  the  out-of-plane 
thruster. 
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Figure  7.12  EPSAT  screen  showing  nozzle  parameters  for  one  nozzle  of  the 
yaw  thruster. 


1 1  9 


A  2  1/2  kg  per  hour  leak  in  one  of  the  thrusters  has  been  reported.  The 
parameters  of  the  1  N  in-line  thruster,  with  a  lower  pressure,  define  a  nozzle 
that  loses  gas  at  this  rate.  The  EPS  AT  screen  describing  this  nozzle  is  shown 
in  figure  7.13.  Note  that  this  nozzle  gives  rise  to  0.5  N  of  thrust.  Figure  7.14 
shows  the  neutral  densities  about  the  subsatellite  due  to  this  leak.  At  the  high 
potential  point,  EPSAT  computes  the  ionization  fraction  to  be  2718,  which  is 
much  greater  than  1 .  Therefore,  the  leak  is  expected  to  cause  a  bulk 
discharge. 
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Figure  7.13  EPSAT  screen  showing  nozzle  parameters  used  to  describe  the 
thruster  leak. 
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